Spin-Orbit induced semiconductor spin guides 
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The tunability of the Rashba spin-orbit coupling allows to 
build semiconductor heterostructures with space modulated 
coupling intensities. We show that a wire-shaped spin-orbit 
modulation in a quantum well can support propagating elec- 
tronic states inside the wire only for a certain spin orientation 
and, therefore, it acts as an effective spin transmission guide 
for this particular spin orientation. 

PACS 73. 21. La, 73.21.-b 



Currently, one of the most challenging issues in con- 
densed matter physics is the injection and control of the 
electronic spin in semiconductor heterostructures. Its 
principal interest comes from the applications point of 
view, since the electron spin in semiconductor nanostruc- 
tures has revealed as a promising candidate to implement 
quantum bits, a necessary ingredient of quantum compu- 
tation jjjj . The feasibility of spin-based electronic devices 
|g|[Q^also relies on the ability to manipulate the spin 
carriers. 

In recent years, spin polarization has been induced in 
semiconductors using optical || and electrical methods 
f|,|. In the case of electrical injection, the technique re- 
lies on external 'spin-alligner' elements, such as magnetic 
semiconductors or ferromagnetic metals, as a previous 
step that polarizes the current which is injected into the 
semiconductor. In this work we propose an alternative 
mechanism that internally selects the propagating spins 
in the semiconductor by means of a space modulation in 
the Rashba spin-orbit (SO) coupling for heterostructures 
with inversion asymmetry. 

We consider a two-dimensional electron gas (2DEG) 
confined to a III-V quantum well whose inversion asym- 
metry produces an electric field in the perpendicular [z) 
direction. As a consequence of the relativistic correc- 
tions this electric field acts on the 2DEG carriers as an 
effective SO coupling field known as the Rashba term 
[@Jl0|. The strength of the Rashba coupling depends on 
the heterostructure's vertical electric field and has been 
shown to be experimentally controllable with a tunable 
gate voltage 11 1 . We take advantage of this tunability to 
define a heterostructure with a space modulated SO cou- 
pling intensity, depicted schematically in the upper part 
of Fig. 1. In practice, the SO coupling variation would 
correspond to a modulated electric field in the vertical 
direction. The well has a constant SO coupling strength 



X R except within a narrow region (guide) of width a 



We model the conduction electrons of this heterostruc- 



ture using the effective mass Hamiltonian 

_ pI+pI , X R (x) 

= 2m* ~1T~ ( Pv<Tx ~ Px<Ty ' 
where the modulated SO coupling reads 
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Equation (Q) also contains the conduction band effective 
mass to* and the Pauli matrices a x and a y corresponding 
to the in-plane electron spin. 

Since the Hamiltonian is translationally invariant in y 
coordinate its eigenstates have well-defined momentum 
in that direction. In this representation the Hamiltonian 
is separable in space coordinates and the eigenstates are 
composed of a propagating longitudinal plane-wave, hav- 
ing y momentum hk, and a spinorial transverse profile 
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where the spinorial part is given in the usual a z basis. 

From the experimental data reported by Nitta et al. 
fblf for an Ino.53Gao.47As/Ino.52Alo.4sAs heterostructure, 
we extract the following parameter values, to be used be- 
low in the numerical applications: to* = 0.05 to c , where 



in. 



A 



where it takes the 'internal' value A^. 



is the free electron mass; X^> = 0.5 x 10~ 9 eVcm; 

R — 1.0 x 10~ 9 eVcm; relative dielectric constant of 
InGaAs e = 13.9; y-coordinate wavevector k = 3.5 x 10 6 
cm -1 , near the Fermi wavevector of a 2DEG with density 
n s ~ 2 x 10 12 cm" 2 . 

Before presenting numerical results with the above 
Hamiltonian and in order to clarify our purpose we con- 
sider the following simplification on Eq. ([!]): let us as- 
sume that the guide can have localized states in the x 
coordinate having a characteristic transverse wavevector 
k x ~ ""/a <C k. Therefore, in a qualitative analysis, 
we can ignore the term containing p x . The eigenstates 
of the resulting approximated Hamiltonian have well de- 
fined spin in the x direction and feel an effective spa- 
tial potential that depends on the <j x eigenvalue, i.e., it 
changes sign from +x to —x spin orientation. Figure 1 
depicts a schematic representation of the up and down 
potentials. It has to be noted that the depth of the ef- 
fective potential well for +x and the height of the barrier 
for —x spin orientation depend both on the spin-orbit 



modulation step AA^ = A^ — A^ ; and the longitudinal 
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y wavevector of the state k. As k increases the attractive 
or repulsive character will also increase linearly. 

The above analysis indicates that transverse space- 
confined modes can exist only for one of the two spin 
orientations thus showing the feasibility of selective spin 
guiding in this kind of heterostructures. We proof below 
that this conclusion remains valid when the Px& v term in 
Eq. (Q) is taken into account. The y-translational invari- 
ance of the eigenstates makes the initial two-dimensional 
Hamiltonian in coordinate space reduce to an effective 
one-dimensional one, corresponding to the transverse 
modes 

h 2 ( d 2 \ ( d \ 

Htr = ~2itf[a*- k ) + Xr{x) + k(Jx ) • 

(4) 

We shall numerically obtain the relevant transverse 
eigenmodes of the Hamiltonian Ti. tr from the resolution 
of the time-dependent Schrodinger equation uniformly 
discretized in the x coordinate and in time. The pro- 
cedure is as follows: we evolve in time an initial spino- 
rial wavepacket and, taking advantage of the harmonic 
time evolution for eigenmodes, we extract the transverse 
eigenspinors and eigenenergies using Fourier analysis of 
the time signals. For a certain value of y-wavevector k 
the time evolution of any wavepacket is decomposed as: 

f ^f^)=^JM^)e-^. (5) 

Besides the eigenenergies cu n k of the different trans- 
verse modes n, Fourier transform at each grid point yields 
the local value of the eigenspinors. This method is appro- 
priate to the present problem since an initial wavepacket 
inside the guide (see Fig. 1) will excite the confined trans- 
verse modes, if they exist, or it will quickly spread to the 
bulk surrounding the guide. We use Gaussian-shaped 
wavepackets in coordinate space with spin oriented in x 
direction, i.e., 

i P (x,t = 0)=^=ef£ X ( ± ) , (6) 
with the <j x eigenspinors 

*£>-( ± l)- CD 

The spatial spread of the wavepacket is chosen close to 
the width of the spin-orbit guide because we want a 
maximum overlap with the confined transverse modes. 
Gaussian wavepackets are appropiate to excite confined 
modes having even parity in the x coordinate while for 
odd parity states it is convenient to use an antisymmetric 
wavepacket. This can be done by multiplying Eq. (|6|) by 
x. 



Physically, the input wavepacket represents an elec- 
tron having well-defined p y and spin oriented along x, 
injected from a lead with transversal dimensions match- 
ing those of the guide. Figure 2 shows the time evolution 
of the density for initial wavepackets with two different 
initial spin orientations: +x (up) and —x (down), using 
a guide width a = 60 nm. It can be seen that for down 
initialization the wavepacket spreads over the transver- 
sal dimension in t ~ 1 ps, with a complete depletion 
of density inside the spin guide. For spin up a part of 
the probability density remains confined to the guide as 
time evolves. This means that our initial up state was 
composed of confined and travelling transversal modes, 
while the down state is only composed of the latter ones. 
For a general orientation of the initial spin the confined 
fraction of the density depends on {<j x ), vanishing if the 
initial spin points in the — x direction and reaching a 
maximum for the +x one. 

To identify and extract the confined modes present 
in the time evolution simulation we use the mentioned 
Fourier analysis technique, focussing on a frequency re- 
gion of the order AX^k. This analysis reveals that for 
our parameter set there is only one confined mode with 
spin polarization mainly oriented in +x direction, as was 
expected from the above discussion. To observe a second 
confined mode (antisymmetric), maintaining the width of 
60 nm, the longitudinal wavevector should be increased 
up to k ~ 6 x 10 6 cm -1 . Actually, for a fixed AA/j the 
number of confined states depends only on the product 
ka. The probability and spin densities of the confined 
mode, related to the wavefunction by 

(p)(x) = \(/) n k-\(x)\ 2 + \(j)nkl(x)\ 2 
(<?z)(x) = \ <Pnk]{x)\ 2 ~ \ Kki{x)\ 2 

(a x )(x) =2Re{<f> nkr (x)4>* nkl (x)} , (8) 

are shown in Fig. 3. 

Probability and spin densities clearly show the confin- 
ing character of the guided mode, the main spin polariza- 
tion being concentrated in +x direction. The mode has 
zero (a y )(x) and a small (a z )(x). This confined mode is 
not an eigenstate of a x , as the transversal Hamiltonian 
Titr does not commute with this operator but, neverthe- 
less, it has an expectation value (a x ) — 0.99, very close 
to the value 1 for eigenstates. 

It is worth to mention that since time-reversal symme- 
try is conserved by the Hamiltonian of Eq. (Q) , Kramers 
degeneracy must hold. Therefore, a counterpart to the 
confined mode discussed above exists. This complemen- 
tary state is also confined to the guide but its spin po- 
larization and longitudinal momentum (j/ y ) are inverted, 
i.e., it is mainly oriented in —x direction and has k' = —k. 

The ability of the guide to confine additional modes 
is enhanced by increasing the guide-width a. Figure 4 
shows the evolution with a, maintaining k to the same 
value as before. Confined modes should lie between the 
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energy edges of the guide, marked in Fig. 4 by dotted 
lines. We note that the second and third modes appear 
for guide widths of a ~ 80 and a w 150 nm, respectively, 
with probability densities having an additional node for 
each successive mode. As the first n — state these 
higher modes have spin almost completely oriented in 
+x direction. 

The confining energies of the guide are of order AA^A;, 
which takes a value of ~ 2 meV for our parameter set. 
Consequently, the practical operation of the proposed 
guide is limited to low temperatures, since the coupling 
with phonons in the meV range would reduce the effi- 
ciency of the guiding effect. In this sense, it has to be 
pointed out that the assignment of a lower spin-orbit 
constant to the guide region is not arbitrary, because, 
as reported by Nazarov et al. spin-orbit coupling 

induces an admixture of pure up and down spin orienta- 
tions in the eigenstates that causes spin decay through 
phonon emission, shown to be the dominant spin deco- 
herence source for quantum dots. To minimize the effects 
of this mechanism the confining structure needs to have 
the minimum spin-orbit coupling. 

In summary, we have shown that a guide structure 
in a semiconductor quantum well, defined by a spa- 
tially modulated Rashba spin-orbit coupling, gives rise to 
transverse-confined and longitudinal-propagating modes 
only for a given spin orientation. This feature charac- 
terizes the proposed structures as spin guides that con- 
stitute natural paths to distinguish and drive the spin 
of the carriers within semiconductors. The feasibility of 
this spin-guiding mechanism is limited to low tempera- 
tures as the energies of the Rashba effect are in the range 
of a few meV's. 
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FIG. 1. Representation of the spatial modulation of 
spin-orbit coupling intensities in the heterostructure contain- 
ing the spin guide. The gray-scale plot shows an upper 
view of the xy plane with the gray tones indicating different 
Rashba coupling intensities (lighter gray corresponds to the 
spin guide). The two lower plots display the effective potential 
in the simplified Hamiltonian, i.e., when p x a y is neglected, for 
+x (up arrow) and —x (down arrow) spin orientations. 
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FIG. 2. Time evolution of the density for an initial Gaus- 
sian wavepacket with well-defined spin orientation for times 
up to ?s 1 ps. Upper panel corresponds to +x initial spin 
direction while the lower one corresponds to — x. The dotted 
lines indicate the position of the spin guide. 
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FIG. 3. Charge and spin densities corresponding to the 
guide-confined mode. As in Fig. 2 the dotted lines mark the 
spin-guide edges. Note the enlarged scale for (a z )(x). 
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FIG. 4. Confined mode energies as a function of the 
guide-width a. Label n indicates successive modes for a given 
a. The dotted lines show the confinement edge energies. Right 
small plots display the probability densities (p) (x) for two se- 
lected modes labeled 'a' and 'b'. 
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